ARKIV FOR FYSIK Band 17 nr 15 


Communicated 13 January 1960 by O. Krery and I. WaLLEeR 


On the foundations of the linear theory of 
irreversible processes. I 


By U. Ustiorn 


1. Introduction 


The purpose of this and a subsequent paper is a study of the linear pheno- 
menological theory of irreversible processes with emphasis on its logical structure, 
its possible probabilistic interpretation, and its relation to statistical mechanics, 
that is to microscopic dynamics. We shall restrict ourselves to systems which 
can be characterized by a finite number of thermodynamic coordinates, the 
kinetics of the system being determined by a set of linear phenomenological 
equations. The thermodynamic coordinates are not assumed to have the same 
parity with respect to a reversion of the direction of the microscopic motion 
and time reversal is introduced into the theory as a linear involutory operation 
on the thermodynamic coordinates. 

Particular interest is devoted to the question of irreversibility. It is pointed 
out that for an arbitrary system of the type considered any initial thermodynamic 
state may be considered as a linear superposition of two states one giving rise 
to strictly irreversible processes and the other to strictly reversible processes. 
The decomposition is unique and very closely related to the possible probabilistic 
interpretation of the theory. It will turn out that reversible processes corre- 
spond to random processes which are of the so called deterministic type whereas 
irreversible processes correspond to random processes which are classified as non 
deterministic or purely random. This connection which is also strongly related 
to the question of macroscopic causality will be discussed in the second part 
of the present work. 

For open systems, under the influence of external constraints which in general 
are time dependent, it is possible to introduce dissipation functions which depend 
only on the state of the system. These functions decrease monotonically in the 
course of time and characterize the approach of the system to a state (in general 
time dependent) which is determined by the external constraints and is inde- 
pendent of the initial state of the system. If in particular the external con- 
straints are time independent, the final state of the system will thus be sta- 
tionary. For systems which obey the principle of detailed balance, that is which 
are described by thermodynamic coordinates which are invariant under time 
reversal, a possible dissipation function is the entropy production. This fact is 
well known as Prigogine’s theorem. On the other hand neither detailed balance 
nor the more general time reversal symmetry is essential for the existence of 


dissipation functions. 
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2. Outline of the phenomenological linear theory of irreversible processes 


Let us consider a closed, adiabatically insulated system under fixed external 
conditions. We shall be concerned only with states of the system which are char- 
acterized by the values of a finite number of extensive variables A = (AP neg AO 
We are thinking of a description of the system in terms of thermodynamic states, 
as outlined in [Uhlhorn, 1960 a]. The statement Att)=a=(e,.%), «") means 
that at the time ¢, the state of the system can be identified with a possible 
thermodynamic state corresponding to a particular constraint chosen from a given 
class of constraints. The particular constraint is determined by the value « and 
conversely the class of constraints determines the possible values of A (t). The 
validity of such a description will not be discussed here. 

There is no further restriction in the assumption that the prescribed constraint, 
that is the statement that the system be closed, adiabatically insulated ete., is 
expressed in terms of the » — m extensive variables Amt). ..., A™(@) in. the form 


A¥ (§) S10! (u=mt+I1, ..., 2) for every time f. (1) 


This can always be achieved by a proper linear transformation of the extensive 
variables. 

We shall assume that under the prescribed constraint (1) the system has ex- 
actly one thermodynamic (equilibrium) state 


Aberg! (iu Lys st), (2) 


The state (2), when once produced, will persist as long as the constraint (1) is 
not changed. This means that if we observe the system, keeping the constraint 
(1) fixed, we will either find that the system ultimately reaches the state (2) 
and no further changes of its state will occur, or no stationary state is reached 
at all. 

It is convenient to normalize the free variables A“ to zero in the equiblibrium 
state (2). Thus in the following we take «“=0 (w=1, ..., m). 

We shall assume that the free variables Al, ..., 4” have a definite behaviour 
under time reversal, or more properly reversal of the direction of motion. We 
shall find it convenient to express this behaviour in terms of a linear trans- 
formation 

To 52 (3) 
which is involutory 
bil Bega oi (4) 


For the given situation, the relevant thermodynamic potential is the, possibly 
generalized, entropy S, which considered as a function S(«) of the free variables 
has a proper maximum in the equilibrium state (2). Defining s,, by 


elaakigng 
she a argenee (5) 


S 1ia the rest of this paragraph the indices ¢, A, 4, v can take the values 1,...,m. Summation: 
convention is assumed for pairs of equal upper and lower indices, 
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=m 
we have the Taylor expansion, valid for small |«|= > ||, 
ual 


S (x) =S (0) —} Sy a a” +0 (| a”). (6) 
The quadratic form s,, ««’ is positive definite: 
a+0O implies s,, a’ a’ >0. (a) 
The entropy function S(«) is invariant under time reversal (3) of its argument, 
and hence 
hp TY Sup = Si,- (8) 
The thermodynamic forces &, corresponding to the thermodynamic coordinates 
a are defined by 
Sy = — Siwy. (9) 
According to (6) 


a os (a) 
oe oil 


+o (|«]). (10) 


We observe that the connection between the variables a” and &é, (9) is rever- 
sible and that both vanish only in the equilibrium state. 

The phenomenological equations, which determine the time development of 
the system are conventionally written in the form (see e.g. [Onsager and Mach- 
lup, 1953]), 


—— al (t) = LM E, (t). (11) 


The phenomenological coefficients are independent of time. They have the fol- 
lowing three properties. The transformation L/”€, is non singular, that is 


I’ —,=0 if and only if €=0. (12) 


This is equivalent to the assumed uniqueness of the equilibrium state. Further 
the quadratic form L/”&,&, is non negative, but not necessarily positive 


L &,€,=3 (LD +L") &, £20. (13) 


This is a consequence of the fact that under the given constraints the entropy 
cannot decrease in the course of time 


d oS da" Ww 2 
se yee Fea b TY é,= TP é,& + >0. 14 
sy esl leap cae ae gu (Meas Carmetl saga (14) 


Conventionally the quadratic expression (13) is called the entropy production in 
the state determined by &. The third property of the phenomenological coeffi- 
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cients is given by the Onsager reciprocal relations 
Ayu _ nw Tua ln = 
Med sits hee Bae, (15) 


We observe that the linear phenomenological theory as formulated above is in- 
variant under arbitrary non singular linear transformations of the thermodynamic 
coordinates. No metric concepts have been introduced so far. However, the 
theory contains two scalar invariants, the quadratic forms corresponding to the 
entropy and to the entropy production respectively. The quadratic form corre- 
sponding to the entropy is non singular and may be used to define a “natural” 
scalar product. Though there is no ambiguity in doing so, we shall find it more 
convenient to use the representation described below. 


3. Representation of the linear theory 


In order to get a clear insight into the logical structure of the linear theory 
of irreversible processes we shall introduce a representation of the theory in 
terms of a space V of “‘state vectors”. Though the physical theory contains 
only real quantities we shall find it mathematically convenient to use a complex 
vector space for its representation. It will turn out that the algebraic structure 
of the linear theory also gives a very strong indication with respect to its pos- 
sible statistical interpretation. 

Let V be an m-dimensional complex vector space of elements wu, v, w, ... with 
a Hermitian scalar product (w|v), linear in its right argument and antilinear in 
its left argument. The norm of a vector v is defined by ||v|]=(v]v)*. Let the 
space V be spanned by m vectors wu’, ..., vw” chosen so that their mutual scalar 
products are real, that is 


(wu | 0”) = (20” | w4).2 (16) 
The reciprocal basis uw, ..., Un is uniquely determined by the relations 
(uM | uy) = df. (17) 


It is an evident consequence of (16) that 
(Un | U,) = (u, | Uy). (18) 
Complex conjugation of an arbitrary vector v= uw" (u,|v) =u, (u“|v) is defined by 
Ov=u" (u,|v) =u" (v|2,). (19) 


According to its definition complex conjugation is an antilinear operation on V 
which leaves the norm invariant. 


1 In this paragraph the indices “4, 7 can take the values 1, ..., m; Summation convention is 
assumed for pairs of equal upper and lower indices. 
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Let the linear operators s, 7, L, M, and Q(t) in V be defined by 


S=| wl) Sy (u"| 


T =|u,) TY (u"| 
L=|u,) L’ (u,| (20) 
Mi Trs 


Q(t)=exp{-tM} (t20). 


We shall call s the entropy operator, 7’ the time reversal operator, and Q (E) 
the semi-group generated by the phenomenological equations (11), for evidently 
Q(t) has the semi-group property 


Q (1) Q (t,) = Q (t, +t) (t, t, 20). (21) 


Further @(t) satisfies the initial condition Q(0)=1. Its infinitesimal generator 
is M=Ls. The operators (20) are real in the sense that they commute with 
the complex conjugation operation C. 

It follows from (5) and (7) that the operator s is positive definite: 


v€V,v+0 implies (v|s|v)>0, (22) 


and hence s', the inverse of s, is well defined. 
The time reversal operator 7’ has according to (4) the involutory property, 
and is related to s by (8). The operator form of these relations is 


P=, (eh =s87, (23) 


where * denotes Hermitian adjoint. 
The Osager reciprocal relations (15) can be expressed in the form 


(FL) = 1 i, (24) 


or equivalently as a condition on the operator M 


M* =(sT) M (sT)", (25) 


that is, M and its Hermitian adjoint M* are connected by a similarity trans- 
formation. 

We shall now turn our interest to the phenomenological equations (11) or 
alternatively to the semi-group of linear transformations determined by Q (t). In 
particular we shall determine the behaviour of Q(t) for large values of ¢. Par- 
ticular attention will be paid to the question whether the time reversal sym- 
metry expressed by (25) represents a real restriction of the mathematical form 
of the phenomenological theory or not. As can be expected the time reversal 
symmetry will not simplify the mathematical treatment except for the particular 
cases in which 7 can be represented by 1 or —1, that is if the variables are 
all of the same kind with respect to reversion of the direction of motion. 
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Let us first recollect the properties we have assumed M to possess: 


a. Non negative entropy production:! s'M*+Ms*>0 
(equivalently M*s+s M>0) 
b. Reality: CM=MC (26) 
c. Onsager’s relations: M*=(s7) M (sT)! 
d. Non singularity: Mv=0 implies v=0. 


To begin with we shall assume only the property a.. This property can also be 
given the form 


Q(t)" sQH=s 0) <s. (27) 


We recall that according to the definition (21) the operator @(t) is a solution 
of the phenomenological equations, with the initial condition Q@(0)=1. The in- 
equality is an immediate consequence of the fact that s(0)=s and that the 
derivative of s(t) is 


a. 
dt 


(t)= —Q (t)* [M* s+sM]Q(t)<0. (28) 


The meaning of (27) is of course that for an arbitrary vector v in V 
Ils? Q (t) o|P=(e|Q (t)* 8 Q(t) |v) < (v| s| v) S| 83 oll’, (29) 


which is a statement of the fact that for an arbitrary initial state v the pheno- 
menological equations (11) cause no decrease of the second order entropy. It 
follows further from (28) that the decrease of s (¢) with increasing t is monotonic. 
On the other hand s(t) is a non negative operator and hence the limit 


s(cc)= lim s (é) (30) 


too 


is a well defined non negative operator. It is evident that s (co) is invariant 
under Q (t): 


Q (t)" 8 (09) Q(t)=s (ce), (31) 
or equivalently 
M* s(co)+s(~)M=0 (32) 
It follows from (29) that 
lim ||s* @ (é) o|| =|]s (co) of. (33) 


Hence we find that 


lim @ (#)v=0 if and only if s(co)v=0, 
t+o 


! An inequality of the form A>B relating two linear operators in V has the meaning that 
A — B is non negative, that is for any v€ V (v| A—B|v) >0. 
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that is, the system can reach equilibrium when started in an initial state v if 
and only if the vector v is an eigenvector of s (co) corresponding to the eigen- 
value 0. 


Let the distinct eigenvalues of M be ju, ..., 4, and the corresponding multi- 
j=k 


plicities m,, .... m, (> mj=m). It is wellknown that the vector space V can be 


decomposed into a direct sum of k subspaces V,, ..., V;. of respective dimen- 
SionS ™m,, ..., m,, Such that each V; is invariant under M and such that (M — 5) 
is nilpotent on JV;.1 The subspaces are in general not mutually orthogonal. 
Taking the initial state v in V; we have identically 


Q (t) v= exp {—t uj} exp {—t[(M—yp,)} o= exp {-ty,}a,([M—p]t)v. (34) 


Here 2;(-) denotes a polonomial of degree less than or equal to m;. From (34) 
we obtain easily some conditions which are necessary for the existence of the 
limit (30). We omit the proofs. 


(a) Re (u;) 20. 


(b) Re (u;)>0 if and only if for every vE€V;,s(co)v=0 or if and 
only if for every v€V;, lim Q(t) v=0. 


too 


(c) If Re (u;)=0 and »v, we V;: (35) 
(v| 8 (t)| w) = (v| s (00) |w) = (v| s| w). 
(d) If Re (u;)= Re (uj) =0, wiApy and vEV;, we V;: 


(v|s (t)|w) = (v| 8 (02 o)| w) = (v|s|w)=0. 


Let W, denote the direct sum of the subspaces V; with Re (y;)>0 and W, the 
direct sum of the subspaces V; with Re (u;)=0. Then V is the direct sum of 
W, and W, and these subspaces are both invariant under M. The unique de- 
composition of an arbitrary vector w€V 


W=Wy + Wy, (Wo E Wi: Wy E W,) (36) 


is in general not a decomposition into orthogonal components. Thus the corre- 
sponding projections Py and P, defined by 


Pyw=W, Prw=w, (37) 


are in general not Hermitian. The invariance of Wy and W, under M can be 
expressed in terms of the projections Py and P;: 


MP,=P,MP,, MP;=P,MP,, (38) 


or equivalently 
Q (t) Po = Po Q (t) Po: Q (t) Pr= PQ (t) Pr (39) 


1 See e.g. [Halmos, 1958 section 58.] 
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The semi-group Q(t) can thus be written as a direct sum of two semi-groups 
Q(t) =Q (t) Po + Q(t) Pr=Qo (ft) + Qi (€), Qo (ty) Or (te) = Or (te) @o (4s) = 9. (40) 

The processes determined by Q, (¢) are strictly irreversible 


lim Q, (t) =0 (41) 


t>00 
and give rise to maximal entropy production 


Q, (0)* 8 Q5 (0) = Po s Py, lim Qy (t)* 5 Qo (t) = 0. (42) 


The processes determined by @,(t) on the other hand are strictly reversible, and 
give rise to no entropy production 


Q7(t)” s Qr(t =PTsP,=P;) 8 (=) P;, (43) 


and wee see that on W, the semi-group Q(t) is equivalent within a similarity 
transformation to a group of unitary transformations 


8? Q,(t)s;* with s;=P;sP,, s7* si =P. (44) 


The operators (44), forming a group of unitary transformations in W;, have 
a spectral resolution of the form 


st Q, (t)s7* => exp {11' t} Ey, 
“ 


where the summation goes over the eigenvalues of the infinitesimal generator 
of (44), that is over the eigenvalues of MP;. The operators H; are mutually 
orthogonal projections forming a resolution of the identity in W,, that is of Py. 
Hence they have the properties 


LEAL =X") 
0 


E* P.=P* E,. 
a" +4)’ A Vig Py E, (45) 


Hence we obtain a spectral resolution of Q;(t) of the form 


Q;(t)= 2 exp {i A’ t} s;* H, st, (46) 


and here s;* HE, s? is the projection into the subspaces V; corresponding to the 


eigenvalue uj;=td' of M in the direct decomposition of V determined by that 
operator. 


We note that the condition (26 a) on M is transformed into itself if M is 
replaced by M* and s by s'. Thus the semi-group Q(t)* is related to d=s7} 
in the same way as @(t) is related to s. In particular we find that 


d(t)=Q (t)d Q(t) (46) 
decreases monotonically with increasing ¢ to a limit 
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lim d(t)=d (co) (47) 


too 


which is a non negative operator bounded by d. The eigenvalues of M* are 
the complex conjugate of the eigenvalues of M and have the same multipli- 
cities. It is thus evident that 


Rank s(co)= Rank d (co). (48) 


It follows that s(co) and d(co) vanish together and in that case Q(t) can deter- 
mine only strictly irreversible processes. If on the other hand the rank of the 
operators is equal to m, the number of dimensions of V, we have necessarily 
s(co)=s and d(cc)=d=s"", and the processes determined by Q(t) are strictly 
reversible. We shall see later that these properties of the semi-group Q(t)* are 
of great value for the statistical interpretation of the linear theory of irreversible 
processes. 

So far we have only taken into account the property (26a) of M. We shall 
now consider the other conditions. The effect of the reality condition (26 b) is 
of course the invariance of the spectrum of M under complex conjugation. 
Further, if V; is the subspace of V which corresponds to the eigenvalue py; 
of M, the subspace CV; corresponds to the eigenvalue &;. M and M”* have the 
same spectrum. 

According to Onsager’s relations (26c) M and its Hermitian adjoint are 
related by a similarity transformation. Thus if V; is the subspace of V corre- 
sponding to the eigenvalue mw; of M, (s7’) V; is the subspace corresponding to 
the eigenvalue yw; of M*. Comparing (27) and (46) we find that s(t) and d(t) 
are related by a non singular transformation 


(sT') d(t) (sT) =s (t). (49) 


Strict reversibility of all processes determined by Q(t) is equivalent to the 


condition 
M*s+sM=0. (50) 


Combining this with Onsager’s relations (26c) we find 
MT =—M, (51) 


Denoting by v(t, %)) = Q(t) v% that integral of the phenomenological equations 
which corresponds to the initial state v) at £=0, we obtain from (51) the relation 


T v(t, 9) =v(—t, Tv). (52) 


The effect of 7 applied to a strictly reversible process is thus really a reversion 


of the direction of motion. 

If on the other hand we assume strict reversibility (50) and the reality con- 
dition (26b) it is always possible to define an operator 7’ which is involutory 
and has the property (51). To see this we note that according to (50) the 
operator s? Ms? is anti Hermitian. Further it is real, that is it commutes with 


C. Hence its spectral resolution is of the form 
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st Ms?= > iAP, (53) 
7 


ae : : = 
with orthogonal projections P, corresponding to the eigenvalues A. With Pa 
also P_, occurs in the expansion (53). The corresponding subspaces have the 
same number of dimensions and hence it is always possible to find a transfor- 


mation 7 which is involutory and maps P,V on P_;V: 
T=), “TP =P (54) 
From (54) and (53) it follows that 


T st Ms} T= —s*t Ms}. (55) 


We may take 7=s } 7s} in order to satisfy (50). 

We shall now discuss the question whether is possible to introduce an ope- 
rator 7 satisfying (26 c) and being involutory if only the conditions (26 a and b) 
are assumed to be satisfied. To facilitate the discussion we shall assume that 
the representative space V has been chosen so that the entropy operator is 
represented by the unit operator. There is obviously no restriction in this choice 
of representation. We note that in this representation a possible time reversal 
operator 7’ will have the properties 


Po Te TOT (56) 


Further L=M, and if this operator is decomposed into its Hermitian and its 
anti Hermitian parts 
L=Ly+ Ly, (57) 


the Onsager relations take the form 

Both Ly, and Ly, are real, that is they commute with C. The first relation in 
(58) shows that 7’ and Ly commute. Hence we may choose a basis in V which 
diagonalizes them both. From (56) it is seen that the possible eigenvalues of 7’ 
are 1 and —1. The condition (26a) is equivalent to the statement that the 
eigenvalues of LZ, are non negative, and imposes no restriction on Ly. The pos- 
sibility of introducing a time reversal operator (56) is thus determined by the 


form of Ly. Let |1), ...,|m) be an orthonormal basis in V which diagonalizes 
Ly. According to condition (26 b) the matrix elements (u|L4|v) of Ly are real 


(u|La|)= — (| Lalu) = (| La]. (59) 


The only possible form of a time reversal operator 7’ is 


P= 2 le)tu(ul  (te=1). (60) 
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In order to satisfy the condition (58) we must have 
Tints (u|L,|v)= —(u|L,4|r). (61) 


It can be seen immediately from (61) that it is always possible to find an 
operator 7 with the required properties if the dimensionality m of V is smaller 
than 3. But for m>3 this is in general not possible and the relation (58) isa 
real restriction of the form of L. 

Finally the condition (26d) is equivalent to the statement that 0 is not an 
eigenvalue of M. This concludes the discussion of the conditions (26) which we 
have imposed on M. 


4. Open system. Dissipation functions 


In paragraph 2 we have considered a closed adiabatically insulated system 
under fixed external conditions. If the phenomenological equations which deter- 
mine the time development of the state of the system are of the strictly irre- 
versible type the system approaches asymptotically the equilibrium state deter- 
mined by the external conditions. This approach is characterized by the mono- 
tonic variation of a dissipation function, which is identical with the first non 
vanishing term in the Taylor expansion of the entropy function in powers of 
the thermodynamic coordinates. We shall now consider the same system under 
conditions which are such that some of the thermodynamic forces have pre- 
scribed values 


£()=&,(),-.4()=4() (<IS<m). (62) 


It is assumed that the state of the system can still be characterized by the m 
thermodynamic coordinates (2) and that the time development is determined 
by the equations 


Fatt” ()= IF, (t = 1 vc HUT, 25M) (63) 
for the “free” thermodynamic coordinates.! In general the constraints (62) will 
prevent the system from being closed and adiabatically insulated. If in partic- 
ular the constraints (62) are independent of the time the system may eventu- 
ally reach a steady state. It is easily seen that in this case the approach to 
the steady state may be characterized by the monotonic behaviour of a dis- 
sipation function, analogous to the entropy function for the closed adiabatically 
insulated system. 

For the discussion of the equations (62), (63) it is convenient to use the 
representation described in paragraph 3. We introduce the, in auto nonortho- 
gonal, projections P, and P,, in V defined in terms of the basis uw’, ...,u™ 


1 In the following primed indices can take the values l,..., 1, doubly primed indices the 
values 1+1,...,m and unprimed indices the values I, ...,m. 
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Ppa =a Do =e e)) (64) 
Pu" =0 (a =U alee) 
P= P,: 
Writing a(t)=a(t)hu,, Et)=E,(t) u", (65) 
we may express (62) and (63) in the form 


; d 
PLE®=-E0), Php aW=PhLE. (66) 


In order to obtain an equation for the time development of the state we use 
the identity 


Peet P. E@)=—{Pie P,Q +Pa oO}, (67) 


to eliminate P,,é(t) from the second relation in (66). We note the existence 
of a unique operator s,, satisfying 


sy Ps Py=Py (68) 

8,, PT =0. 
It follows from (68) that s,, is non negative. Hence (67) has the unique solution 
Py E(t)= — sy [Phe P, E+ Pha @)=—sy [st EQ +a(t)], (69) 


and we obtain from (66) the equation 


d ; ' 
ap Ph, a(t)= —Ph,Ls,, a(t) + Ph Dé (t)—s,, s+ € (t)). (70) 
Introducing 
M,,=PLLS M7 )= Peli @— see (71) 


we obtain the equation (70) in the form 
d 
a [Pi «())= —M,, (Pi a«(t)+ Pi 7 (t)). (72) 


From the non negativity of D+ L* follows that of 
Pe (Tits Ve ate LPaget Pmt oP alee 
and hence 87 Uo Mine = 0, (73) 


That is, we find that M,, and s,, are related in the same way as M and s. 
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Both s and s,, are non negative operators, but S;, 18 In general singular. It is 
now evident that 


5, (4) =exp {-t My} $1, @Xp { it M,;} (74) 


decreases monotonically from the value s,, at t=0 to some non negative limit 
81,(c°). The properties of s,,(co) are determined by the spectrum of pee dtn 
particular we note that if MW can give rise only to processes of the strictly ir- 
reversible type we have s,,(cc)=0. 
The formal integral of the equation (72) is 


t 
Phi a(t), = exp {— (§—t,) M,,} Ph «(t,) — i exp {—(t-t’) M,,} M,, Pint’) dt’, (75) 


to 


the first term on the right hand side representing the changes which are due 
to processes inside the system and the second term representing the changes 
due to the external constraints. If in particular s,, (co) =0, that is if the changes 
taking place inside the system are strictly irreversible, and if 7 (t) is sufficiently 
wellbehaved as a function of ¢, for instance if ||7 (t)||<C@< co for all values of 
t, we can let f)+— oo in equation (75) and obtain the “long time” solution 
of (72) 


Pia(t)e=— | exp {-t M,,} UW, Pin t-t')dt (76) 


0 


that is the solution which any solution of (72) approaches asymptotically. We 
want to characterize the approach of an arbitrary solution «(t), determined by 
e.g. «(0) to the “long time” solution, which is independent of «(0), by the 
monotonic decrease of a dissipation function. The difference between the two 
solutions (75) and (76) is 

0 


PY; [a (t)g — & (t) oo] = exp {-£M} gy | exp {t’ M,,} My, Ph n(t')dt'] (77) 


and 
0 


|| si, PTs [oe (9 — & (t) oo] ||? = || S41 (2)? [x (0) — J exp {t’ M,,} M,, Pin (t')dt’] 


—o 


? (78) 


decreases monotonically to zero. Hence (78) which evidently is a unique func- 
tion of the state «(t) of the system, may be used to characterize the approach 
towards the “long time’ behaviour, which is only dependent on the external 
constraints. This dependence is causal in the sense that the state at the time 
t is determined by the values of the external parameters at times prior to ¢. If 
in particular the external constraints are time independent, 7 (t)=%, (78) can 


be written 

Il st: [a (o— Pr nll? (79) 
and the steady state is given by 

Pi Asteady —Ph, Y- 
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The existence of the dissipation function (79) may be considered, at least par- 
tially, as an affirmative answer to a wellknown question raised by the Ehren- 
fests [Ehrenfest, 1911 p. 82] about the possibility of characterizing the steady 
state by the minimum of some function which could be considered as a gener- 
alization of Boltzmann’s H-function. 

It is interesting to compare the characterization of the steady state discussed 
above with Prigogine’s theorem of minimum entropy production [Prigogine, 
1947]. To obtain Prigogine’s theorem we observe that a steady state, given the 


constraint P,é=€ is characterized by 


d ° 
eect Ostendy = 1a ecto — i Ps Eseenay = Ph LP; €=0. (80) 


Assuming that this equation has at least one solution P,, &steaay, and this will 
always be the case if Z is non singular we may describe an arbitrary state by 


E=Eteaay tA &, Pi NESE. (81) 
Consider now the expression (€|Z|£&) and apply (80). 
(€| L| &) = (Ectoaay | L| Eatenay) + (Af | Z| A €). (82) 


In the case that L is Hermitian (82) represents the entropy production in the 
state determined by &. The quadratic form is non negative and evidently takes 
its minimum value for Aé=0. Thus, if LJ is Hermitian and non singular there 
will always exist a unique steady state determined as the state of minimum 
entropy production. This is the content of Prigogine’s theorem. It is obvious 
that under the conditions of Prigogine’s theorem the system will reach its steady 
state in the course of time. 

To see that the decrease of the entropy production to its stationary value is 
monotonic we write the integral of the equations (66) in terms of the thermo- 
dynamic forces, observing that the constraints are independent of time, 


E@j=exp {—t[s, PLE) Pf O46, Pa OP) 4s, P| 
=(e(Ph DP Ay ey PRE taney ee 8 ee) 
Here (s,, Pi, LP,,)"* is defined by 
6) PREP Heer Py SP (84) 
Pie, Py DP, {)=0. 


Let A be an arbitrary operator in V and consider the derivative of the matrix 
element (A€(t)|A|A E(t) =(AE(t)| Ph AP, |AE(): 


d 
ap SeMl4lAE)= 


—(Aé()|PREP 8, Ph APgHPYAP,, Sk LP A Aid) Bo) 
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and this expression is non positive for every choice of Aé(t) if and only if 
(2b Pa)éa (Pid Pa (Pea Peel Pa Ue, eS 0. (86) 


We note that the solution (83) is valid under the sole assumption that 
8, Pi LP,, is invertible in the sense of (84). Hence it is valid whether L is 
Hermitian or not, whether L is of the strictly irreversible type or not. How- 
ever, €steaay Will be the limit of (83) only if all the processes taking place in- 
side the system are strictly irreversible. 

Assuming now with Prigogine that LZ is Hermitian and taking 4=Z in (85) 
we find that (86) is satisfied, and comparing with (82) we find that the decrease 
of the entropy production is monotonic. 

Another possibility of satisfying (86) is to take 4=s~'. Using the definition 
(68) of s, we find that (86) becomes 


PL IAP PoP Ped DP, 20, (87) 


and this inequality is valid on the basis of the non negativity of D+ L*. Hence 
(AE|s"*| A) is a dissipation function of the type required by Ehrenfest. It is 
easily seen that if in (A&|s *|A&é) the thermodynamic forces are expressed as 
functions of the thermodynamic coordinates we obtain the dissipation function 
(79) for the particular case of time independent constraints. 

It is thus seen that the dissipation function (79) is different from Prigogine’s 
entropy production but that the monotonic decrease of both functions is due 
to the same property of the phenomenological equations. 

The dissipation function (79) has been defined in terms of the operator 5,,, 
given by (68). We shall conclude this section with a discussion of the thermo- 
dynamic signification of this operator. We recall that the entropy operator s is 
defined by 

es 
Oa" Oa” 


s=—|u") (0) (u” |, (88) 


where S(«) is the entropy function. The inverse of s can be expressed in terms 
of the Massieu-Planck function [Uhlhorn, 1960 a] 


S 
¥ (B)=B, a! +S (a) (#.-25 (0), (89) 
52 
Pn ae Ht (90) 


=| Uy) 5B, OB, (0) 


oy 
0 Bur 8 By 


Hence Pee gee) 0) (wy |. (91) 


Expressing s,, in the form : , 
Saas | a!" ) Stay (W (92) 


we find that sy) = S114 is determined by 
eo. 


6.0 (93) 
Su1Ay é By 6 By yh 
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We recall that doubly primed indices take the values +1, ...,m. Considering 
the generalized thermodynamic potential Y, defined by 
Peace See 
ve = Y; (Das ar) Bi ou e809 gett) — S rc S (94) 
pial (alog 


we find [Uhlhorn, 1960 a, equation (53)] that 


a2 
ey, 
Sti — Sars me Bar 


the derivative being taken at the reference state (2). Introducing A a=" (t))— 
— a (t)oo, (77), we find that the dissipation function (78) can be written 


Cie io - 
—_—,—— A Noe 96 
Oa 6a” ‘ i ee 


Hence we have found how the dissipation function is related to the thermo- 
dynamic functions of the system. Returning to Ehrenfest’s question, we compare 
the expression (93) with the generalized form of Boltzmann’s principle [Uhlhorn, 
1960 a, (71)]. We find that for time independent constraints the dissipation 
function (93) determines the probability distribution of the fluctuations in the 
stationary state. 

Dissipation functions of the type discussed above can be defined for con- 
tinuous systems also as far as the kinetic equations can be considered as linear, 
that is for systems in mechanical equilibrium, if the phenomenological coefficients 
are independent of the space coordinates. 
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